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Abstract — Indoor localization systems that are accurate and 
robust with respect to propagation channel conditions are still 
a technical challenge today. In particular, for systems based on 
range measurements from radio signals, non-line-of -sight (NLOS) 
situations can result in large position errors. In this paper, we 
address these issues using measurements in a representative 
indoor environment. Results show that conventional tracking 
schemes using high- and a low-complexity ranging algorithms 
are strongly impaired by NLOS conditions unless a very large 
signal bandwidth is used. Furthermore, we discuss and evaluate 
the performance of multipath-assisted indoor navigation and 
tracking (MINT), that can overcome these impairments by 
making use of multipath propagation. Across a wide range 
of bandwidths, MINT shows superior performance compared 
to conventional schemes, and virtually no degradation in its 
robustness due to NLOS conditions. 

I. Introduction 

Indoor localization systems face many difficulties when they 
are to be deployed in practice (|T|. Optical systems might 
become unusable in emergency situations due to fire and 
smoke; signal strength based systems are sensitive to fading 
effects; fingerprinting systems need dedicated training phases; 
etc. Range-measurement-based systems using radio signals, 
which are in the focus of this paper, offer advantages like the 
capability of being integrated in existing radio devices such 
as smartphones, the potential of low-power implementations 
and moderate infrastructure requirements. However, they face 
challenging performance impairments caused by propagation 
effects like strong reflections or diffuse scattering. 

To counter these impairments, a large signal bandwidth is 
beneficial for two main reasons: First, the temporal resolution 
and therefore the attainable accuracy is increased. Second, the 
line-of-sight (LOS) path can be separated from other signal 
parts easier, which improves the robustness, i.e. the percentage 
of cases in which the attainable accuracy can be achieved. For 
these reasons, systems using ultra-wideband (UWB) signals 
are promising candidates for both accurate and robust indoor 
localization |2|. 

In this paper, we examine the performance of a conventional 
range-based UWB indoor localization and tracking system in a 
representative indoor environment with several anchor nodes. 
Using channel measurements along a reference trajectory, 
mimicking the walking path of a pedestrian user, we show 
the dependence of the robustness and accuracy of the local- 
ization on changing propagation conditions. The LOS path 
is estimated using the maximum-likelihood method as well as 
using a low-complexity ranging algorithm. These estimates are 



provided as inputs to an extended Kalman filter (EKF), which 
tracks the user's position. The presence of strong reflections in 
the measurement signal can lead to a large bias in the estimate 
of the propagation delay of the direct path between user and 
anchor This is especially true in non-LOS (NLOS) situations, 
where the direct path is blocked. Overall, a degradation of the 
robustness of the system is the result. 

Usual countermeasures are to detect such situations based 
on signal statistics and discard the corresponding measure- 
ments f3l, (4). In advanced schemes like f5l, f6l, machine 
learning techniques are used to mitigate the range bias in 
NLOS situations. However, these approaches discard position- 
related information that is inherent in the geometric structure 
of reflected MFCs. The usefulness of this information has been 
shown in llTl-lfTOl. 

We have shown in previous work JQl- lfTTI that by using 
additional floor plan information, these reflected signal paths 
can be used for localization. Extracting these paths from 
the channel measurements and associating them to the floor 
plan, signal reflections can be seen as direct signals coming 
from virtual anchor (VA) nodes. In this way, one physically 
existing anchor node is turned into a set of virtual anchors, 
and NLOS situations no longer render this anchor useless for 
localization. This scheme, which we call multipath-assisted 
indoor navigation and tracking (MINT), can thus deliver the 
desired robustness with respect to the propagation conditions. 
However, the effect of diffuse scattering is especially pro- 
nounced in indoor environments due to the presence of objects 
like, e.g., furniture. This makes the extraction of the reflected 
components from the measurements and the required data 
association challenging. To highlight these issues in a practical 
setting, this paper extends our previous work in [11] with the 
use of multiple base stations, a detailed performance analysis 
of the proposed algorithms and a comparison to standard 
tracking schemes employed in indoor localization. 

The key contributions of this work are: 

« We show the performance and the robustness issues of 
a conventional, EKF based tracking system with respect 
to channel conditions. This is done in a representative 
indoor environment using measured signals. 

• Making explicit use of multipath, robustness and accuracy 
can be increased. We discuss an implementation of the 
proposed MINT approach and examine its performance. 
Results show the excellent robustness of MINT with 
respect to NLOS situations. 

« The results obtained highlight general practical consider- 



6 






1 












1 1 


- 






' ^'' 






4 










' ^72 






^^^-== 




BS3 












' J13 Virtual signal^ 


^ 


2 


=s>^^- 




140 1 




-Traiectorv 

1 
2 






— 











I 


220 180 ^^ 


■^ 


i^' 


■ ( 






















>< 






^"^ 


"-V^ 


J8 ^y^ 


* 

V 
















118 












>. 

■S 












_ 










411 




X 


■> 




















^ 
^ 
«* 






^ 


•>» 


•S. 






-4 




•* 


^ ^ 266 

' K .■■ ■■ 












On, 

*>. . 

in. 


" 


-6 


412^^ 




1 


1 










"l^22 





-5 



10 



15 



20 



25 



x[m] 



Fig. 1. Indoor tracking scenario considered in this paper The black dotted line illustrates the user's trajectory in the scenario, selected position indices are 
given to allow easier inteipretation of later results. Blue circled crosses ai'e the base stations, while black and grey squares with crosses denote some few 
illustrative VAs for BSs 3 and 4. The area shaded in light gray denotes the geometric visibility region of the LOS path of BS 4, illustrating the NLOS region 
caused by the hypothetical obstruction in the center that is used for performance evaluations. Additionally, the thin and bold blue lines show two example 
propagation paths that can be modeled by VAs and that are potentially visible in the NLOS region. 



ations for indoor localization systems, i.e. the importance 
of a large bandwidth and propagation channel conditions 
is shown. 

This paper is organized as follows: Section HIl explains the 
scenario, the measurement campaign and the geometrical mod- 
eling of the environment. A short overview over performance 
bounds for MINT is given in Section |llll while Section |Ty] 
discusses channel estimation and ranging algorithms. Tracking 
is explained in Section [V] and a detailed discussion of results 
is given in Section IVll 



II. Scenario, Measurements, and Geometric 
Modeling 

Fig. [T]shows the scenario used in this paper. It is the ground 
floor of a large three-storey building with open ceilings in 
the corridor areas in all but the uppermost floors. Walls are 
made of reinforced concrete (shown as black outer lines) and 
doors (shown in grey) are made of metal. The long grey outer 
lines on the upper and lower side are large windows with 
some metal pillars in-between. We placed four base stations 
(BSs) at known locations and measured the UWB channel 
between them and a moving agent. The latter was moved along 
a trajectory consisting of 220 points spaced by 10 cm. The 
obstruction shown near the center of the corridor has not been 
present in the measurements. It is introduced artificially in 
the performance evaluations to show the influence of NLOS 
scenarios on the localization, c.f. Section [VT] 



A. Geometric Modeling of the Environment 

If a floor plan of the building is available, the BS positions 
Sii, i — 1, ... ,4, can be mirrored with respect to reflecting 
surfaces like walls to obtain so-called first-order virtual an- 
chors (VAs) at positions a'^. Using these, the procedure can 
be repeated to obtain second-order VAs a"i, and so on. The 
visibility regions of these VAs, i.e. the regions where the 
corresponding reflections are possible, can be computed by 
ray tracing as explained in [121 . If the corresponding MFCs 
are detectable at the moving agent, they can be associated to 
the VAs and used for localization. Fig. [T] shows a subset of 
VAs for BSs 3 and 4 together with two example signal paths. 

In an analogous manner to [111, the overall set of VAs for 
the i-th BS is denoted as 



J^^ = {j 



J 



1, 






(1) 



Here, the superscript indicating the VA order has been dropped 
as it is not important in the upcoming discussion. A function 
/vis(aj , p) determines the visibility of the j-th VA at a position 
p, i.e. 

^ , , I 1, if VA a, is visible at p 
/v.(a„p) = |^^ else. ^'^ 

The set of expected visible VAs at a specific position p^ can 
then be computed evaluating (|2]i for all VAs in ([TJ 



A. 



(») 



{^ 



,a^^(.)} = {aj : /vis(aj, p^) = 1}. (3) 



In the following, we omit the BS index i wherever it should 
be clear from the context that the respective quantity exists 
for all BSs. 

It can be seen from (|3]l that a set of anchor nodes sig- 
nificantly larger than the set of BSs can potentially used for 
localization. We call this approach multipath-assisted indoor 
navigation and tracking (MINT). However, the challenges are 
as manifold as the potential gains: First, the MFCs correspond- 
ing to the VAs need to be detected in the signal received at 
the agent, which depends on influences such as e.g. material 
parameters or reflection angles. Second, together with deter- 
ministic MFCs, also some scattered, non-specular components 
will cause a certain level of clutter in the measurements. Third, 
the detected MFCs need to be associated to the VAs to make 
their position-related information accessible. Summarizing, a 
moving agent needs to perform detection of MFCs in the 
presence of possibly strong diffuse multipath and noise as well 
as subsequently a robust data association procedure to infer its 
location. 

B. Channel Measurements 

Measurements were performed in the frequency domain 
with a vector network analyzer (VNA) over the full FCC band- 
width from 3.1 to 10.6 GHz. These measurements have also 
been used in ITOl - llTSI . The discrete frequency representation 
of the channel's transfer function at position p^ is denoted 
as Hi[k]. Using a frequency spacing of A/, it corresponds 
to a Fourier series (FS) of a time-domain channel impulse 
response (CIR) hi{T). This CIR is periodic with a period of 
T'max ~ 1/^/^ which is the maximum resolvable delay. The FS 
representation allows an arbitrarily fine time resolution Ar to 
be achieved. Using an IFFT with size N-pYi — \{^I^'r)~'^^, 
the time domain equivalent baseband signal is obtained as 

-i27r(/,-/o)t 



r,{t) ^WFlN^{H(,[k]S[k]}^' 



(4) 



In ^, S[k] is the discrete frequency domain representation of 
a raised-cosine pulse s{t) with pulse duration Tp and roll-off 
factor /3r [13|. It is used to cut out the desired band of the 
signal as well as to filter it to facilitate an MFC estimation 
algorithm. The frequencies fc and /o are the desired center 
frequency and the lower band edge of the extracted band, 
respectively. This procedure is similar to [14|. Examples for 
ri{T) obtained in this way are shown in Fig. |2] for BS 3 over 
the trajectory. 

III. Channel Influence on Localization 

The received signal at the i-\h trajectory position p^, 
corresponding to the processing in (3), is modeled as ||9] 

ri{t) = ^ ak,es{t - Tk,i) + s{t) * Vi{t) + wi>{t). (5) 
fc=i 
It consists of a sum of Kg deterministic MFCs, scaled and 
shifted by the complex amplitudes ak,i and propagation delays 
Tfcf, respectively. The signal vi{t) is a stochastic process 
modeling the diffuse multipath (DM), i.e. diffuse scattering at 
rough surfaces and other propagation mechanisms that are not 




Fig. 2. Received signals r^ (t) in logaiithmic scaling over the trajectory for 
BS 3. The white dashed lines indicate expected delay paths of MFCs modeled 
by VAs up to order two, computed by geometric ray-tracing. Indices of some 
of the VAs shown in Fig. ^ are given. 



captured by the specular part of the channel. It is observed at 
the agent as a convolution with the transmit pulse s{t). Finally, 
wt{t) denotes additive white Gaussian noise (AWGN) with 
two-sided power spectral density Nq/2. 

We expect that the deterministic part of rg{t) can be 
modeled by the set of VAs Ae- In f9], we have derived the 
Cramer-Rao lower bound (CRLB) on the position error for the 
MINT scheme, based on the received signal ri{t). It is given as 
the inverse of the equivalent Fisher Information matrix (EFIM) 
fS]. If all Ki deterministic MFCs are orthogonal at p^, i.e. 
there is no path overlap in the delay domain, the EFIM can 
be decomposed as 

o^2a2 Ke 
Jp, = -^ Y. SINRfc,,Jr(<^fe,£) (6) 



where /3 denotes the effective (RMS) bandwidth of s{t) and 
c is the speed of light. The decomposition is instructive since 
it separates the channel's influence from the geometry of the 
VAs. The former is given by the signal-to-interference-plus- 
noise-ratio (SINR) of the deterministic MFCs. For the fc-th 
MFC at position pi, it is defined as ||9l 

SINRm = ^Jri\ V (7) 

Here, interference is understood as the contribution of the 
DM, which is given by the power delay profile (FDF) of the 
DM, Si,,i at the delay Tk^e of the respective MFC, scaled with 
the pulse duration. This FDF is defined as the second central 
moment of the DM process i'i{t). 

The influence of the geometry is accounted for by the 
ranging direction matrix Jr(0fc,^) LSJ, li9J, which is defined 
as 

'''^'^- e) cos{<pk,i) sin{(l)k,i) 



Okl 



cos 



'^k,l) 



cos[<pk,e) sm((^ 



(8) 



This matrix has one eigenvector pointing in the direction <j>k,£, 
the angle from the VA at a/^ to position pi. The SINR in O 
is used in the performance evaluations to emulate changes of 
propagation conditions caused by the artificial obstruction that 
is shown in Fig. [T] 

IV. Channel Estimation and Ranging 

Conventional time-of-anival (ToA) localization requires an 
estimate of the time-of-flight of the signal from the BS 
to the agent. This process is called ranging [15|. In this 
work, we consider two different ranging algorithms, maximum 
likelihood (ML) and jump-back search-forward (JBSF). The 
former is computationally intensive since it requires channel 
estimation, while the latter is a popular low-complexity al- 
gorithm. The MINT approach requires range estimates to all 
the MFCs that are modeled by VAs, i.e., an estimation of the 
deterministic part of the UWB channel as for ML ranging. 

A. Channel Estimation and ML Ranging 

It has been shown that the MFCs corresponding to the 
VAs can carry a large fraction of the energy of the received 
signal 1(121 . hence we extract those Ki MFCs that result in the 
minimum energy of the difference signal 



{h./} = arg min 



nit) 



E 

fe=i 



ak,i>s{t - Tkj) 



dt (9) 



where T is the observation interval. Of course the parameter 
Kg should be related to Kg in (|5]). However, without the 
knowledge of the position, a value has to be chosen that allows 
for the extraction of the geometrically relevant paths over the 
entire scenario. The choice of this parameter will be discussed 
in Section [VLB] 

For this work, we assume a separable channel, i.e. there 
is no overlap between deterministic MFCs. In this case, (|9]l 
can be performed path-per-path. The amplitudes ak^g are 
nuisance parameters to this estimation and can be obtained 
by a projection of the unit-energy pulse s{t), shifted to the 
estimated delay ffe £, onto the received signal 



&k,i = / rl{t)s{t - Tk,i)dt 



(10) 



where * denotes complex conjugation. With the separable 
channel assumption, this procedure is similar to lfT6l . 

For MINT, the set of estimated arrival times in (|9]l, con- 
verted to distances, constitutes the measurement input to the 
tracking algorithms. Additionally, a data association (DA) 
scheme is necessary, which will be described in Section IV-AI 
For ML ranging, only the propagation delay of the first arriving 
path, the LOS path, is interesting and can be obtained from 
dg as 

^Los f = ™in{ffe^£}. (11) 




T[ns] 



Fig. 3. Exemplary ranging results for BS 4 and position I = 120 (c.f. Fig. 
^, using Tp = 0.5 ns. The lower plot is a close-up. The black line denotes 
the received signal, while the red signal illustrates the reconstmction using 
the ML channel estimator The latter also provides estimates for the MPC 
arrival times r^ ; (grey dashed lines), where the ML range estimate (in red) 
is the minimum of these. Here, JBSF (light grey) outperforms ML due to path 
overlap. The simplified implementation of the ML channel estimator assumes 
a separable channel and only detects the larger second component. 



B. Jump-back Search- forward (JBSF) Ranging 

If real-time localization is desired, a channel estimation 
might not be feasible. JBSF is 



like (|9]l might not be feasible. JBSF is a threshold-based 
ranging algorithm that also works for low-complexity receiver 
structures like energy detectors |T5l. Backwards from the 
maximum of the received signal at delay r^ax, a search-back 
window of length Tsb is defined. The JBSF range estimate is 
the minimum delay within this window that exceeds a certain 
threshold 



.-(JBSF) 
^LOS,f 



l{7 



- Tsb < r < Tn, 



\ri{T)\>A^}. 



(12) 
The choice of the threshold A^ is usually based on the noise 
level. 

An illustrative example for the channel estimation and rang- 
ing algorithms is given in Fig. [3] In this case, the LOS path is 
not separable from the rest of the channel. Since the amplitude 
of the second arriving MFC is larger, it is detected in an earlier 
iteration of the ML estimator Due to the separable channel 
assumption, the ML estimator is constrained not to estimate 
any MFCs within one Tp of already detected ones, hence 
the LOS path is not detected. However, the LOS component 
is above the threshold of JBSF and within its search-back 
window, so JBSF can outperform the ML estimator in this 
case. 

V. Tracking and Data Association 

The conventional tracking schemes based on ranging and the 
MINT approach are both implemented in this work using ex- 
tended Kalman filters (EKFs) flTl . flSl employing linearized 
observation equations for the estimated path delays. As in IfTTI . 
we use a simple constant-velocity state-space model of the 



moving agent node. Stacking the position p — [px,Py]'^ and 



velocity v = 
model reads 



Vx,V, 



]-^ in the state vector x — [ 



pT^^TjT jjjg 



Xi+i ==Fxe + Gna 



(13) 



1 AT 
1 AT 
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Xf 
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1 
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AT" 
2 







AT. 



na,^. 



Here, AT is the update interval and ria.^ denotes the driving 
process noise, which allows for motion at non-constant ve- 
locity. The measurement inputs for the conventional tracking 
scheme using EKFs are either the ML or JBSF ranging 
estimates. Using (fTTI ) or ( fT2] i. ranging is performed to each of 
the four BSs. The entries of the vector of distance estimates 



are obtained as 2: 
the BS index. 



(«) 



CT, 



(ML,i) 



LOS J 



or z 



(«) 



CT, 



(JBSF,? 
LOS,£ 



, where i is 



A. MINT Using EKF with Data Association (DA) 

For the MINT scheme implemented using an EKF, the MPC 
estimation in (|9]l only results in an unordered set of distance 



Ai) 



^i) 



estimates, i.e. Z^ — c-{ff, [}, where the cardinality of the set 

is \Zy\ — K^ ' . To obtain a suitable input for the observation 
equations of the EKF, a DA scheme is necessary. As the 
association to the BS is known, we only consider the DA 
for one BS and drop the index i. After the DA, all associated 
distance estimates are stacked into one measurement vector. 

The employed DA approach is the same as used in 1 1 1 1 
and is just shortly re-sketched here. We introduce a set 
Ci of correspondence variables [19], whose n-th entry C£„ 
represents the association of the n-th entry of Zi and is defined 



Ci,n 



I j, if Z£.„ corresponds to VA a^ (j = 1, . . . , Nva) 

10, if Zi n corresponds to clutter 

(14) 

Using the predicted position of the EKF, p^ , the set of 
expected visible VAs Ae can be computed using (|2]i and ([3]). 
From this, the set of expected path lengths, denoted as Vi, can 
be obtained by calculating the Euclidean distances d{a.j,pj) 
between all VAs a^ G Ai and the position pj . The task of the 
DA is now to match the two sets Zg and Vg in an optimum 
way. 

The approach used is optimal sub-pattern assignment [20 1. 
Noting that the cardinality of Vg is Kg (c.f. Q) and assuming 
that Ke > K(, which can always be ensured by filling up 
Zi with dummy clutter, we search the sub-pattern of Zi that 
fits best to the set Vg. Using a distance function d*^'^'^'(-), that 
simply saturates at the cut-off distance dc, this is expressed as 

Kt 

TTopt = arg min Y] d^'^'=^((i£,i, z^^^J. (15) 

Here, 11^ is the set of permutations of positive integers 
up to Ki and the vector TTopt holds the permutation with 
the minimum cumulative distance in the sense of the cutoff 




Fig. 4. 
for all 



200 220 



Instantaneous K-factor with respect to the LOS component -R'los 
BSs over the trajectory using a pulse duration Tp = 0.5 ns. 



distance function d^'^''>{-). Hence, the first Ki entries of TTopt 
contain the indices of those entries of Zi that have been 
optimally assigned to 2?£. 

As a last step in the DA process, we reject those measure- 
ments for which d^'^''>{di^i,Zi^Tr„p,i) — dc, i.e. those that have 
been assigned at a distance greater or equal than the cutoff 
distance. This is done to have a mechanism to tune the DA to 
the uncertainty in the MPC delay estimates. As UWB signals 
are used, if an MPC is detectable at all, its delay will have 
a relatively small uncertainty that depends on the used pulse 
duration Tp and the inaccuracies of the floor plan. Hence, the 
correspondence variables in (IT4l are obtained as 



Q,»i — 



if vTj = n and d'-'^'^^idi^, z^^^^^^J < 4 
else. 



(16) 



The optimal sub-pattern assignment can be efficiently imple- 
mented using the Munkres algorithm 121]. 

B. Genie-aided Data Association (GADA) for Evaluation 

The DA discussed above depends crucially on the EKF's 
predicted position pj. To enable an evaluation of the local- 
ization performance where the errors due to false DA are 
mostly eliminated, we use a so-called genie-aided version of 
the DA (GADA). The GADA reUes on the true position p^ for 
performing the DA, i.e. in the calculation of the set Vg. This 
can be used to benchmark the usefulness of the MPC delay 
estimates for locaUzation using MINT. 

VI. Results 

This section illustrates the accuracy and robustness of the 
proposed indoor tracking algorithms for the scenario in Fig. [T] 
A first rough evaluation of the propagation conditions can be 
obtained by evaluating the K-factor with respect to the LOS 
component, i.e. the ratio of its energy to the energy of the rest 
of the CIR. It is estimated using (O and ( ITOl ). where only the 
LOS component is subtracted from ri{t), i.e. 



K, 



"LOS,. 



LOS.^ 



/y \ri{t) - aLOS/s(^ - Ti^os,eWdt 
where we note that the pulse s{t) is energy-normalized. 



(17) 



TABLE I 

Parameters and results for the performance evaluation. For the performance measures, a pair of values is given where the first 
one indicates the scenario without the obstruction and the second one with the obstruction. 





Tp = 0.2 ns 


Tp = 0.5 ns 


Tp = 1 ns 


Tp = 2 ns 


Tp = 4ns 


Process noise variance (t„ [m^/s*] 


""a = ("zAT )^' '^o'variance matrix Q = cr^GG"^, with Vx,m-a = 1.5 m/s 


EKF measurement noise variance a^ [m^] 


0.01 


0.01 


0.04 


0.04 


0.09 


DA cutoff distance dc (MINT) [m] 


0.3 


0.3 


0.5 


0.5 


0.6 


JBSF SB-window [ns] 


100 


JBSF tlireshold i [1] 


0.4 


0.4 


0.3 


0.3 


0.3 


Mean number of assoc. MPCs (EKF-GADA) [1] 


18.2 / 17.6 


18.6/ 18.2 


18.3/ 17.9 


15.0 / 13.9 


11.9/ 10.2 


Mean number of assoc. MPCs (EKF-DA) [1] 


18.1 / 17.5 


18.5 / 18.2 


18.2/ 17.7 


15.1 / 14.0 


11.8/ 10.1 


RMS position en'or MINT (EKF-GADA) [m] 


0.03 / 0.032 


0.039 / 0.041 


0.06 / 0.069 


0.076 / 0.083 


0.117/0.136 


RMS position error MINT (EKF-DA) [m] 


0.073 / 0.087 


0.067 / 0.065 


0.118 / 0.197 


0.122 / 0.149 


0.207 / 0.304 


RMS position error EKF, ML ranging [m] 


0.05 / 0.052 


0.058 / 0.065 


0.079/0.113 


0.129 / 0.149 


0.353 / 0.396 


RMS position error EKF, JBSF ranging [m] 


0.113/0.15 


0.351 / 0.418 


0.243 / 0.297 


0.249 / 0.34 


0.406 / 0.714 


average HDOP MINT (EKF-GADA) [1] 


0.41 /0.39 


0.41 / 0.38 


0.39 / 0.39 


0.43 / 0.43 


0.47/0.51 


average HDOP MINT (EKF-DA) [1] 


0.67/0.73 


0.56 / 0.53 


0.56 / 0.65 


0.58 / 0.61 


0.68 / 0.76 


average HDOP EKF ML ranging [1] 


1.19/1.21 


1.23/ 1.2 


1.35/ 1.23 


1.09 / 0.96 


1.04 / 0.8 


average HDOP EKF JBSF ranging [1] 


1.23/1.28 


1.27/ 1.52 


1.56 / 1.49 


1.28/1.24 


1.03/ 1.14 



Fig. |4] shows i^LOS.f for all BSs over the trajectory, es- 
timated using the original measurements (without the ob- 
struction) with Tp ~ 0.5 ns. Some geometric influences, like 
the drop in iiTLos in the NLOS regions of BSs 3 and 4 
or the general distance dependence can be easily observed. 
Other causes of large variations are not immediately intuitive 
and require more careful inspection of the local propagation 
conditions. For example, the drop of K^los for BS 1 around 
position 90 is due to strong scattered components from metal 
pillars within the lower windows, c.f. Fig. [T] 

A. Evaluation Setup 

The "obstruction" introduced in Fig. [T] is meant to model 
a change in the propagation environment, e.g. a group of 
people. As it causes NLOS regions for all BSs, its influence 
on the tracking performance is interesting. MINT is not 
aware of this obstruction in its floor plan knowledge. The 
original channel measurements where performed without this 
obstruction, hence they need to be modified. For this, we 
calculate the positions which are now in NLOS due to the 
changed conditions. At these positions, we reduce Ku:>5 by 
a factor of 10 dB, which is done by subtracting the LOS 
component and adding it again with a suitable amplitude and 
its original phase. To enable a fair comparison with MINT, 
we also change the amplitudes of all the MPCs affected by 
the change in visibility conditions. We do so by reducing their 
SINR (Q also by 10 dB. Note that this factor leaves the LOS 
component still detectable in most cases, while for the later 
arriving MPCs it causes an effective loss of the respective 
component 

B. Choice of Parameters 

a) Measurement parameters: We perform the perfor- 
mance evaluation for several exemplary signal bandwidths, i.e. 



the pulse duration of the raised cosine pulse is chosen as Tp S 
{0.2ns, 0.5ns, Ins, 2ns, 4ns}. The corresponding — 3dB 
bandwidths are {5 GHz, 2 GHz, 1 GHz, 0.5 GHz, 0.25 GHz}. 
Using a roll-off factor of /3r — 0.5, the center frequency is 
always fc = 7 GHz, except for Tp = 0.2 ns which corresponds 
to the full measurement bandwidth, hence the overall center 
frequency fc = 6.85 GHz has to be used. Table |I] contains an 
overview of all relevant parameters for these setups. 

b) EKF parameters: The EKF's measurement noise co- 
variance matrix is R = cr^I, where I is the identity matrix. 
For the variance ct^, we choose values based on the expected 
ranging uncertainty, i.e. it is slightly increased with the pulse 
duration Tp. The same consideration applies for the DA cutoff 
distance dc used in (fTsT i and ( fTSI l. Table U lists the values used 
in the different measurement parameter setups. 

The choice of the process noise variance cr^ depends on the 
maximum velocity in e.g. the s-direction, denoted as w^.max- 
The process noise, assumed to be Gaussian, should allow 
movement at this velocity. Hence, we select the 3ct point of 
the velocity noise as Vx.nax and obtain for the corresponding 
process noise variance in the acceleration domain 

In this paper, we choose to model pedestrian motion and set 
Wo;. max — 1.5 m/s. The update interval is AT = 0.1 s, resulting 
in a magnitude of the velocity in x- or y-direction of 1 m/s 
along straight segments of the trajectory. The conventional 
tracking schemes and MINT always use the same a1 and a^ 
for every pulse duration considered. 

c) Initialization: All EKFs are initialized with the correct 
position and zero velocity. We note that this is a more chal- 
lenging problem for MINT than for the conventional schemes, 
since the DA depends on the estimated position. However, we 
have addressed this problem for a different implementation 
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Fig. 5. CDFs of the ranging en'or for ML (red) and JBSF (blaclc) with and 
without the obstruction for all BSs and Tr, = 1 ns. 



of MINT in f22] using a Gaussian-sum filter Il23l . In this 
approach, a bank of parallel weighted Kalman filters (KFs) 
is used to represent multiple initial position hypotheses. After 
convergence, the KFs with small weights can be discarded. 
This can also be applied for the MINT implementation in this 
paper 

d) Channel estimation and ranging: The choice of K( in 
^, i.e. the number of extracted components, is based on the 
average number of expected visible VAs in the localization 
environment. For the performance evaluations, we take into 
account VAs up to order two and fix K£ = K — 20. 
However, we define a threshold between the noise floor and 
the maximum amplitude of the received signal, below which 
no further MFCs are extracted. The noise signal ■wi{t) is 
estimated from the pre-LOS part of rg{t). Denoting the time 
average of a signal with (■), the threshold is obtained as [12] 

A^^e = 7 • (max{|r,(t)|} - {\Mt)\)) + {\Mt)\)- (19) 

In this paper, 7 = 0.1 is used. 

The JBSF ranging threshold A^ in (fTZt is also chosen using 
a relative threshold ^, in the same manner as 7 in (fT9] l. For 
every pulse duration used, we selected a threshold at which 
ranging works well in the unobstructed scenario, c.f. Table IT] 
The length of the search-back window Tsb is chosen for each 



BS as the maximum distance between the BS and any first- 
order VA, expressed in ns. This corresponds to the maximum 
delay difference of the LOS component and any first-order 
reflection. The reasoning behind this is the assumption that 
if the LOS component is not the maximum itself, then the 
maximum will be among the first-order MFCs. A value of 
100 ns is selected for the geometry in Fig. [1] It approximately 
corresponds to the excess delay of the reflection path on the 
far right end of the corridor 

C. Ranging Performance Results 

A comparison of the performance of ML and JBSF ranging 
algorithms in terms of the ranging error CDF is shown in 
Fig. |5] with and without the obstruction. A pulse duration 
of Tp = 1 ns has been used for this. The most striking 
result is that JBSF closely approaches ML ranging for all 
cases, at a much lower computational complexity. One reason 
for the good performance of JBSF is the rather low noise 
level in the VNA measurements. However, when looking 
at the instantaneous ranging errors over the trajectory (not 
shown), some isolated outliers of JBSF can be observed. The 
implications for the localization performance will be discussed 
in the next subsection. 

D. Localization and Tracking Performance Results 

Table |I] contains detailed performance results in terms of 
RMS position error, the mean number of associated MFCs 
for MINT and the average horizontal dilution of precision 
(HDOP). The latter value is calculated as the average of the 
instantaneous HDOP values over the trajectory. These are the 
ratio of instantaneous position error and RMS ranging error 
IfTSl, hence they provide snap-shots of the quality of the 
localization geometry. For MINT, the RMS ranging error is 
computed from those range estimates that have been associated 
to VAs by the DA procedure. 

Looking at the RMS position error in Table U obtained using 
both ranging methods reveals that the EKF using ML ranging 
still clearly outperforms the one using JBSF. This behavior 
would not be suspected by looking at the ranging CDFs alone. 
The few isolated outliers of JBSF mentioned above cause the 
gap in localization performance. As this subsection shows, 
MINT does not show such robustness issues with respect to 
ranging outages. 

Fig. |6] contains the position error CDFs for all considered 
pulse durations. MINT using EKF-GADA always shows the 
best performance of all tracking algorithms. Although this 
is just an upper bound for this implementation of MINT, it 
shows the potential gain of exploiting multipath propagation. 
This is not only true for the achieved accuracy, but also for 
the robustness with respect to difficult propagation conditions, 
since MINT using EKF-GADA shows almost no performance 
loss for the obstructed scenario. Also with realistic DA, MINT 
is able to overcome impairments caused by the obstruction 
better than the conventional approaches. The conventional 
approaches show significant performance impairments, espe- 
cially when the pulse duration is increased. 
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Fig. 6. CDFs of the position error over the trajectory for all considered pulse 
durations Tp. The left column shows runs without the obstruction (c.f. Fig.[T). 
in the right column, signals have been modified according to the obstruction. 



Only for pulse durations of Tp ~ 0.2ns and Tp = Ins, 
MINT with DA can not provide a performance gain. In the first 
case, the channel estimation can resolve very closely spaced 
MFCs that are challenging for the DA due to their clustered 
structure. In the second case, interference by pulse overlap 
situations causes the detection of several MFCs to fail. Since 



neither the channel estimation (|9]l nor the DA as proposed here 
can deal with path overlap situations, this is identified as an 
important topic for future research. 

An unexpected result is obtained for Tp — 4 ns, which, 
corresponding to a — 3dB bandwidth of 250 MHz, is no longer 
UWB by definition [24 1. Although MINT generally benefits 
from a large bandwidth, it can yield the largest performance 
and robustness gain here. The conventional tracking schemes 
show outliers due to very large ranging biases, while MINT 
still exploits information from some MFCs, although their 
number decreases significantly (c.f. Table |I]i. The mentioned 
large ranging biases are also the reason for the lower HDOF in 
the obstructed scenario in these parameter setups. Furthermore, 
it should be noted that in terms of the RMS position error, 
MINT using the EKF-DA always outperforms the conventional 
EKF approach using JBSF ranging. 

An illustrative example for Tp = 0.5 ns and the NLOS 
scenario is shown in Figs. |7]and[8] The estimated trajectories 
of the agent node are displayed in Fig. Q While MINT and the 
EKF using ML range estimates can both follow the trajectory 
well, the EKF employing JBSF ranging is severely impaired 
by the NLOS regions of BSs 1 and 4 (at the start of the 
trajectory), BS 4 (around position 120, c.f. Fig. [1} and BS 3 
and 2 (near the end of the trajectory). As Fig. |8] confirms, in 
the first two of these NLOS regions, rather large values of the 
instantaneous HDOF occur This highlights the bad geometry 
for conventional schemes when there are ranging outages. For 
MINT, most of the HDOF values along the trajectory are below 
one, indicating excellent geometry for localization. This is due 
to the large set of virtual BSs spanned by the VAs (c.f. (|6) 
and (EJ). 

VII. Conclusion and Future Work 

We have presented a comparison of conventional indoor 
UWB ranging and tracking schemes and an implementation 
of the proposed MINT approach using data from a measure- 
ment campaign in a representative indoor environment. The 
influence of NLOS situations on the accuracy and robustness 
of these systems has been investigated. Results confirm the 
potentials of MINT: Using floor plan information, multipath 
propagation can be efficiently used to obtain position-relevant 
information in situations where conventional approaches fail. 
Even at comparably low bandwidths, MINT shows a clear gain 
in localization performance. 

Ongoing and future work includes the modeling of floor 
plan uncertainties and the combination of localization and 
channel estimation. The latter can be used to reduce the 
amount of clutter in the extracted path delay estimates. Also, 
other implementations of MINT, e.g. using soft classification 
in the data association, are within the scope of our work. 
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